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We present the first remotely operated vehicle investigation of megabenthic communities (1004-1695 m 
water depth) on the Hebrides Terrace Seamount (Northeast Adantic). Conductivity- temperature-depth 
casts showed rapid light attenuation below the summit and an oceanographic regime on the flanks 
consistent with an internal tide, and high short-term variability in water temperature, salinity, light 
attenuation, aragonite and oxygen down to 1500 m deep. Minor changes in species composition (3-14%) 
were explained by changes in depth, substratum and oceanographic stability, whereas environmental 
variability explained substantially more variation in species richness (40-56%). Two peaks in species 
richness occurred, the first at 1300-1400 m where cooler Wyville Thomson Overflow Water (WTOW) 
mixes with subtropical gyre waters and the second at 1500-1600 m where WTOW mixes with subpolar 
mode waters. Our results suggest that internal tides, substrate heterogeneity and oceanographic interfaces 
may enhance biological diversity on this and adjacent seamounts in the Rockall Trough. 



The ecology of seamount communities in the deep sea is governed by multiple drivers including envir- 
onmental heterogeneity, intrinsic community dynamics (migration, extinction, species interactions 1 ' 2 ), 
and disturbance (fishing) history 3 , the importance of which varies between seamounts. The presence of 
internal tides, heterogeneity in seabed substrata, and temperature changes often vary with depth on a single 
seamount, creating substantial changes in species turnover 4,5 and genetic variability 6 even on the same seamount. 

Although such general trends in the ecology of seamount communities in the deep sea are emerging, continued 
seamount exploration and monitoring have global importance in an era of expanding human activities penetrat- 
ing deeper into the oceans. Despite a long history of scientific exploration in areas such as the Rockall Trough in 
the northeast Atlantic 7 , only 0.4-4.0% of global seamounts have been explored with scientific purpose 8 . For 
example, biological communities on only two of the three seamounts in the Rockall Trough (Rosemary Bank and 
Anton Dohrn Seamounts) have actually been surveyed. This hinders our understanding of how these seamounts 
contribute to the wider structure and function of deep-sea ecosystems in the northeast Atlantic 8 . Unexplored 
seamounts are also at risk of having an under-developed evidence base to underpin their protection because the 
biological or ecological value of these seamounts, necessary for conservation prioritisation 9 , cannot be 
ascertained. 

Today, dropframe camera and remotely operated vehicle (ROV) technologies can provide contiguous and geo- 
referenced replicates of high-resolution imagery ideally suited to seamount exploration and examining paradigms 
of seamount ecology across space and time in the deep sea. Dropframe camera surveys of the Anton Dohrn and 
Rosemary Bank Seamounts in the Rockall Trough have revealed a range of habitats including cold-water corals 
reefs, coral gardens, and dense aggregations of habitat-forming sponges and xenophyophores 1013 , but these also 
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showed evidence of fisheries impacts 14 . Benthic communities on the 
third seamount in the Rockall Trough, The Hebrides Terrace 
Seamount (HTS: Figure 1), were only recently surveyed using a cam- 
era-equipped ROV combined with oceanographic profiling 15 . The 
HTS is an elliptical flat-topped guyot formed in the early Cenozoic 
during three phases of volcanic activity beginning 62 Myr after the 
mass extinction event at the end of the Cretaceous 16 . The HTS is 
surrounded by Plio-Pleistocene debris deposits up to 900 m thick, 
but is also marked by evidence of more recent seismic events along its 
southern edge 17 . Numerous canyons, gullies and escarpments occur 
along the southern and western flanks of the seamount 18 . It rises from 
depths of around 2300 m to 1000 m deep 19 , several hundred metres 
deeper than the summits of the other seamounts in the Rockall 
Trough. 

The HTS was already known to be important for both marine 
cetaceans and fish, with recent sightings of fin whale Balaenoptera 
physalus, long-finned pilot whale Globicephela melas, harbour 
porpoise Phocoena phocoena and beaked whale, likely Mesoplo 
don bidens 20 . Large spawning aggregations of orange roughy 
(Hoplostethus atlanticus) on the HTS were discovered in the early 
1990s but these were rapidly overexploited within a decade and the 
fishery has since been abandoned 21,22 . Dense occurrences of cold- 
water corals and sponges at the other two seamounts, combined with 
occurrences of marine cetaceans, fish spawning aggregations, and 
possibly an isolated stock of H. atlanticus on the HTS, suggest that 
oceanographically dynamic conditions at this seamount 23 may help 
concentrate high quality food sources for both visiting and resident 
seamount fauna 24 , a functional role that should also support biologic- 
ally rich seamount communities. 

The oceanographic setting of the Rockall Trough bathing its sea- 
mounts is also both spatially and temporally dynamic 25,26 , with the 
interface between water masses being important zones of biotic 
change 7,27 . Surface waters (SW) are derived from the westward- 
flowing North Atlantic Current (NAC) and overlay an upper layer 
centred around 600-1200 m of northward-flowing warm saline 
Eastern North Atlantic Water (ENAW) that forms in the Bay of 
Biscay. In the vicinity of the HTS, ENAW is often mixed with 
mid-depth outflow of Mediterranean Water (MW) and Sub- Arctic 
Intermediate Water (SAIW) the latter which derives from the west- 
ern boundary current of the subpolar gyre, a mixture that varies 
annually 26 . Below 1200 m, cooler dense overflow waters are derived 
from the Norwegian Seas that flow through the Faroese Channels 
and over the Wyville Thomson Ridge 28 . This Wyville Thomson 
Ridge Overflow Water (WTOW) is recorded as far south as 55°N 
in the Rockall Trough and is itself mixed with MW and SAIW, 



although WTOW signals >1500 m water depth are likely attributed 
to large periodic overflow events 28 . Beginning around 1500-1600 m, 
WTOW transitions to Labrador Seawater (LSW) down to about 
2000 m water depth 26 . Cool, fresh LSW is formed by deep winter 
convective mixing in the Labrador Basin and exported eastwards by 
the subpolar gyre into the southern entrance of the Rockall Trough 2 '. 

Thus, the objective of our study was to use the ROV images and 
oceanographic profiles to investigate environmental variability, spe- 
cies richness and megabenthic species composition across the HTS. 

Results 

Oceanographic and hydrographic variability at the HTS. ROV- 
based conductivity-temperature-depth (CTD) measurements showed 
that the HTS is bathed by several water masses (Figures 2 and 3, 
Tables 1 and 2). The summit at (—1000 m water depth) is mostly 
covered by ENAW, but between 800-1000 m there are distinct 
waves of high salinity MW. ENAW extends to about 1200 m, 
where mixing with WTOW and SAIW occurs. The transition 
from ENAW to WTOW at about 1200 m overlaps with high 
variability in the particle attenuation coefficient (Cp), salinity and 
temperature (Figures 2 and 3, Table 2) that persists until about 
1500 m. Between 1200-1600 m the waters were mostly WTOW, 
which is cooler and fresher, less dense and with higher dissolved 
oxygen (DO) concentrations than ENAW (Figures 2 and 3). 
Deeper than 1500 m, there is less variation in temperature, salinity 
and DO (Table 2) and thus we term this layer as oceanographically 
more "stable" than the "dynamic" layer above 1500 m water depth. 
Near the end of the dive transects at depths s 1600 m there was 
increasing penetration of LSW characterised by even higher DO 
(>260 iimol/L) than WTOW (Figures 2 and 3; Table 2). 

DO, dissolved inorganic carbon (DIC), Cp, total alkalinity (TA), 
aragonite saturation state (f2) and pH all varied with depth (Figure 2, 
Table 2). Steep declines in the concentration of oxygen occurred at 
800-1000 m water depth with an oxygen minimum of 203 umol L -1 
at 1000 m water depth, after which DO steadily increased up to 
256 umol L _1 at 1930 m water depth. Oxygen saturation was nearly 
100% at the surface, decreased to about 68% near 1000 m and 
increased below 1000 m back up to about 77% at 1930 m (Table 2). 

The seamount summit experienced little change in Cp, ranging 
between 0.15-0.16 irT 1 (Figure 2). Thereafter sustained high vari- 
ability in Cp ranging between 0.1-0.2 m _1 persisted down to 1500 m 
water depth (Figure 2, Table 2). 

Aragonite saturation decreased from 1.33-1.38 at the seamount 
summit to as low as 1.11 at a depth of 1930 m (Figure 2, Table 2). 
Estimated O. aragon it e for the ROV transects revealed saturation states 
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Figure 1 | Location of the Hebrides Terrace Seamount in the southcentral Rockall Trough and the three ROV dive transects on the seamount summit 
and flanks (inset). Maps were generated using ArcGIS version 9.3. 
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Figure 2 | Depth profiles from the full CTD cast (0-1930 m water depth), showing measured (A): temperature (°C), (B): salinity (psu), (C): dissolved 
oxygen concentration (DO), (D): particle attenuation coefficient (Cp), (E): dissolved inorganic carbon (DIC), (F): total alkalinity (TA), and 
calculated (G): aragonite saturation state, and (H): pH (total scale). In figures (E-H) black crosses indicate the mean value of the data collected from water 
samples (error bars are 1 standard deviation). 



< 1 on several occasions, reaching as low as 0.76. Q ara gomte > 1 is a 
state when thermodynamically, aragonite mineral precipitation 
occurs, whereas Q < 1 is a state when thermodynamically aragonite 
dissolves. Thus our CTD measurements indicate a degree of periodic 
undersaturation on the deep seamount flanks. 

Natural history. A total of 109 megafaunal taxa were identified 
(Supplementary Table SI). The giant agglutinated xenophyophore 
Syringammina fragilissima was the most frequently observed species, 
present in 66% of images. Elasmobranch egg cases from at least two 
species were observed at 16 stations on the flanks in waters 1523- 
1651 m deep; some were observed among coral branches of the reef- 
framework-forming scleractinian Solenosmilia variabilis, similar to 
the way in which shark egg capsules are observed on cold-water coral 
reefs formed by another scleractinian, Lophelia pertusa 30 . Their 
identities have not yet been confirmed with archived museum 
collections so these egg cases were not included in the benthic 
analyses. One sponge species identified tentatively as 7 .Halichon- 
dria (Eumastia) sitiens from dive 037 is known to contain bioac- 
tive anti-cancer compounds 31 . 

A total of 26 taxa were tentatively identified to species level, all of 
which have broad North Atlantic or global ranges. All four fish 
taxa identified to species level are distributed across the North 
Atlantic including the oreostomatid Neocyttus helgae that frequently 



co-occurs with seamount cold-water corals and sponges 32 . Many had 
Lusitanian or Arctic affinities that suggest the HTS is inhabited by a 
rich mixture of faunae from the north and south via the different 
water masses that bathe this seamount. We caution however that 
because our species identifications were based solely on ROV images, 
some species could actually represent cryptic species. Morphological 
and molecular data as presented in Cardenas et al 33 for the geodiid 
species in the Atlantic boreo-arctic region will be needed to validate 
our identifications. 

Species richness. In contrast to low species numbers on the summit 
(dive 036, mean number of species per transect = 4.5, s.e.m. = 0.1), 
the seamount flanks (dives 035 and 037) were densely colonised by 
species-rich assemblages (mean = 8.7, s.e.m. = 0.3 and mean = 7.1, 
s.e.m. = 0.2, respectively) including habitat-forming fauna such as 
Syringammina fragilissima, geodiid sponges and the reef framework- 
forming scleractinian coral Solenosmilia variabilis. Megafauna inha- 
biting Solenosmilia framework on the HTS were similar to those 
living among Lophelia framework on nearby banks and seamounts 
e.g. black corals of the genus Stichopathes, glass sponges, encrusting 
sponges, ophiuroids and crinoids 11-3436 . 

Although there was a lot of variability in the number of species 
from image to image, richness on the flanks appeared to begin 
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Figure 3 | Water mass properties at the HTS. (A): T-S and DO conditions 
encountered on the HTS through a vertical profile from the surface to 
1930 m. (B): T-S and DO conditions for the ROV dives 035, 036 and 037 
with density (a) contours. Water mass indicators for SW, (surface waters), 
ENAW (Eastern North Atlantic Water), WTOW (Wyville Thomson Ridge 
Overflow Water), MW (Mediterranean Waters), SAIW (Sub-Arctic 
Intermediate Water) and LSW (Labrador Seawater) properties are shown 
as either boxes or mixing-lines. Representative values for water masses 
were taken from McGrath et al. 26 for the Rockall Trough. 

increasing at around 1250 m water depth, peaking around 1300- 
1400 m, then decreasing again (Figure 4). A second peak in richness 
seemed to occur around 1500-1600 m (Figure 4). Non-parametric 
regression of species richness on the flanks showed that a substantial 
proportion of variation in the number of species (between 40-56%) 
was explained by depth, substratum, water mass identity, and inter- 
actions between them (goodness of fit R 2 dive035 = 0-55, R 2 dive037 = 
0.40). Little variation in richness on the summit was explained by this 
combination of variables (goodness of fit R 2 dive036 = 0.02). 

Changes in species composition. Substrata and oceanography did 
not vary at the summit as they did on the flanks, the latter which 
exhibited significant changes in species composition related to 



stronger downslope environmental gradients (Figure 5). Strong 
community dissimilarity was detected from an analysis of simi- 
larities (ANOSIM) between communities on the summit and the 
flanks, with habitat- forming sponges, corals and xenophyophores 
down the flanks but which were virtually absent on the summit. 

Community dissimilarity on the northeast flank (dive 035) was 
strongly related to downslope shifts in the predominant substrata, 
which changed from coarse sand to gravels, cobble, bedrock and 
coral framework (Supplementary Table S2; Figure 5). Communi- 
ties significantly (ANOSIM global R = 0.32, p = 0.001) differed 
across most substrata except for cobble and very coarse gravel, and 
cobble and bedrock, each of which could not be significantly discri- 
minated (p = 0.40 and 0.89, respectively; Supplementary Table S2). 
With increasing depth, similarity of percentages (SIMPER) analyses 
revealed how assemblages changed from one dominated mostly by 
sponges such as 7 .Periphragella lusitanica, Hamacantha sp. and 
Hexadella dedritifera on very coarse gravel, to S.fragilissima, sponges 
and ophiuroids on cobble, then to one characterised by sponges, 
ophiuroids and the antipatharian coral Stichopathes cf abyssicola 
on bedrock, and coral framework characterised by Solenosmilia var- 
iabilis, sponges and ophiuroids (Supplementary Table S2). 

Community dissimilarity on this flank was also strongly related to 
differences in oceanographic variability (ANOSIM Global R = 0.17, 
p = 0.001; Supplementary Table S2). SIMPER detected differences 
between assemblages bathed by more dynamic versus stable waters 
including more occurrences of sponges, antipatharians, solitary cor- 
als, and the reef-building coral S. variabilis at depths > 1500 m and 
more Syringammina fragilissima at oceanographically more 
dynamic stations (Supplementary Table S2). 

Community dissimilarity on the seamount's southeastern flank 
(dive 037) was also related to downslope changes in substrata but 
not as strongly as for dive 035 (ANOSIM Global R = 0.15, p = 0.001; 
Supplementary Table S2). Shallow stations began on coarse gravel 
fades and changed to patches of cobble, bedrock, and coral frame- 
work then finally to fine sands at the deepest stations (Figure 5; 
Supplementary Table S2). Communities that did significantly differ 
were characterised by SIMPER as being comprised of sponges, 
Stichopathes cf abyssicola and Solenosmilia variabilis on coarser 
grained substrata to those with higher occurrences of Syringam- 
mina fragilissima on finer grained sands (Supplementary Table S2). 

Community dissimilarity on dive 037 was also related to down- 
slope shifts in oceanographic variability (ANOSIM Global R = 0.14, 
p = 0.001; Supplementary Table S2). Differences identified by 
SIMPER were driven primarily by higher proportions of Sticho- 
pathes cf abyssicola and ophiuroids at deeper more stable stations 
while the large glass sponge IPeriphragella lusitanica occurred more 
frequently at oceanographically dynamic settings <1500 m deep 
(Supplementary Table S2). 

Although community dissimilarity was closely related to envir- 
onmental variability, in fact only a modest amount of change in 
species composition could be directly explained by gradients in these 
variables or by the effects of geographic distance between stations. In 
total, environmental and geographic distance explained statistically 
significant (p < 0.01) amounts of change in species composition on 
the seamount flanks (12-14% in dives 035 and 037, respectively) and 
on the summit (3% in dive 036; Supplementary Table S3). The 
change in substrata and oceanography during dives 035 and 037 
detected during the SIMPER analyses were again reinforced in the 
species-environment biplots as were the species associations includ- 
ing habitat-forming corals characterising deeper stations (Figure 6). 

Discussion 

Our study established for the first time the high benthic biological 
diversity and occurrences of habitat-forming megafauna on the HTS 
alongside an oceanographically and hydrographically dynamic set- 
ting. Changes in species composition were related to differences 
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Table 1 ROV transects on the Hebrides Terrace Seamount together with ROVCTD range for temperature, 
and estimated aragonite saturation state 

Dive 035 Dive 036 


salinity, and estimated oxygen % 
Dive 037 


Date of dive (dd/mm/yyyy) 


09/06/2012 


10/06/2012 


1 0 /OA/901 9 

1 \J / VJO/ 1 Z 


Transect length 


3157 m 


1040 m 


9079 m 


Subset # of images analysed 


319 


84 


9Q9 


Subset duration (hours:minutes:seconds) 


06:56:30 


01:42:57 


OA-9A-9/1 
UO.ZO.z4 


Geomorphology 


flank 


summit 


rlank 


Depth range 


1238-1695 m 


1004-1021 m 


1234-1661 m 


Temperature range 


4.56-6.40 


7.50-7.62 


3.97-5.43 


Salinity range 


34.97-35.17 


35.22-35.24 


34.95-35.09 


Oxygen saturation (%) 


63.6-83.8 


58.3-74.3 


62.3-82.0 


Aragonite saturation state 


0.77-1.40 


0.87-1.36 


0.77-1 .35 


Total # observed megafauna species 


80 


33 


56 



between the summit and flank, which differed in the presence of an 
internal tide impinging on the sloping flanks of the seamount at 
1200-1600 m water depth. Oceanographic variability and changes 
in substratum were related to changes in species composition on the 
flanks, but neither variable directly explained many of these changes. 
In contrast, between 40-56% of changes in species richness was 
explained by this environmental variability. Notably, two peaks in 
species richness centered at the interface between different water 
masses suggest a critical role of oceanographic variability in enhan- 
cing richness. Thus our findings reinforce the notion that seamount 
fauna do respond to environmental variability across different spatial 
scales, e.g. summit versus flank or downslope, but in species-specific 
ways 37 and in such a way that they cannot be directly explained using 
categorical variables such as water mass, or predominant substratum. 

Community dissimilarity between the seamount summit and 
flanks was likely related to seamount-generated flow variability such 
as internal tides and mixing 1,38 . Critically, internal tides and en- 
hanced mixing bring food particles to deep-sea benthos that 
increases community similarity across wide areas over which these 
features occur 39 . Hydrographic variability creates community dis- 
similarity on the Great Meteor Seamount off the Azores 40 , and our 
evidence also supports this notion for the HTS where long sections of 
the flanks were inhabited by large habitat-forming sessile suspen- 
sion-feeding corals and sponges. We note that very little evidence of 
fishing disturbance was evident during the dive on the summit: thus 
it is more likely that community dissimilarity between the summit 
and flanks were likely due to the presence of hydrographic features 
such as an internal tide. 

Surface currents are enhanced and strongly baroclinic over the 
HTS 41 , the abrupt topography of which generates an internal tide 41-42 
that is strongly displaced in two zones centered on the 600 m and 
1200-1600 m water depth contours 43 . Given the high density and 
diversity of large suspension-feeding fauna strictly on the flanks, our 
biological and CTD data support the notion of this internal tide 
delivers higher concentrations of particles to the seamount flanks 
(Figures 2, 3) than at the summit, which in turns substantially alters 



species composition between summit and flanks. Benthic megafauna 
on the adjacent Anton Dohrn Seamount are remarkably similar to 
those documented here, with enhanced densities of habitat-forming 
sponges, corals and xenophyophores at comparable depths 13,44 . 
Notably, there is also a strong internal tide that propagates from 
Anton Dohrn 41 , which likely plays an equivalent role on this sea- 
mount as it does on the HTS. 

Community dissimilarity was also detected between oceanogra- 
phically dynamic versus more stable water mass layers. Our CTD 
measurements documented oceanographic transitions from warmer 
more saline ENAW to deeper layers of WTOW and then to cooler 
and fresher LSW. This transition structures benthos across the wider 
Rockall Trough 7,27 as similar interfaces do for seamount taxa glob- 
ally 4,45 . Of particular note at the HTS were the relatively higher fre- 
quency of Solenosmilia variabilis at depths with more short-term 
thermohaline and oxic stability, and more frequent occurrences of 
the large xenophyophore Syringammina fragilissima at shallower 
stations that had higher short-term oceanographic variability. 

Some of the species-specific responses to downslope changes on 
the HTS can be attributed to changes in chemical oceanography. For 
example, the oceanographic setting at depths beyond 1500 m water 
depth is ideally suited to the environmental niche of S. variabilis. 
This species can inhabit cooler less saline waters that are less satu- 
rated in aragonite and higher in oxygen 46 than another ubiquitous 
reef-forming coral in the Rockall Trough, Lophelia pertusa, which 
tends to inhabit seamounts < 1200 m water depth in the Rockall 
Trough 44 . Other variables such as DIC, alkalinity, pH and potential 
water mass density constrain the distribution of L. pertusa* 7 , however 
effects on occurrences of S. variabilis are not unknown. Depth is 
known to be a strong predictor of suitable habitat for S. variabilis 
as is the presence of a seamount 46,48 . Given that S. variabilis also 
exists along restricted salinity and temperature ranges 46 , our study 
suggests that the relationship with increasing depth on seamounts is 
perhaps better explained by oceanographic stability and geochemical 
constraints such as oxygen concentration, which may also be depth- 
related. 



Table 2 | Means and ranges of oceanographic variables on the HTS recorded by the full CTD cast. Ranges of salinity, temperature, particle 
attenuation coefficient (Cp), oxygen and aragonite saturation states were two to three times higher at stations 1200-1500 m (the 
oceanographically "dynamic" layer) throughout most of the WTOW layer than at (oceanographically stable) > 1 500 m water depth 





Salinity (psu) 


Temperature (°C) 


C p (m-') 


O2 saturation (%) 


Aragonite saturation 


Means 












1200-1500 m 


35.07 


5.38 


0.14 


74.56 


1.18 


1500-1700 m 


34.97 


4.27 


0.16 


77.80 


1.18 


Ranges (max-min) 












1200-1500 m 


0.08 


0.93 


0.06 


3.63 


0.04 


1500-1700 m 


0.03 


0.40 


0.03 


1.35 


0.03 
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Figure 4 | Depth-related trends in species richness per image during dives 035 and 037. Blue bars indicate what appear to be peaks in species richness. 



On the continental slope to the west of Britain, the agglutinated 
xenophyophore Syringammina fragilissima is typically encountered 
in depositional environments characterised by enhanced surface 
production and particle flux 49 . Thus enhanced frequencies of this 
species at oceanographically variable stations on the flanks of the 
HTS are probably explained by the greater range of Cp values 
recorded by CTD at depths 1000-1500 m which would supply more 
particles to these stations. 

Substratum heterogeneity significantly altered species composi- 
tion on the flanks of the HTS as it does on other seamounts 1,4 . 
Occurrences of deep-sea reef framework- forming corals in particular 
alter species composition and effectively triple species richness 50 . 
Coral reef framework on the HTS was formed by Solenosmilia var- 
iabilis, and comprised a very rich mixture of associated species 
including black corals, bamboo corals, cup corals, sponges, crinoids 
and tunicates, all of which further enhance diversity themselves by 
providing more habitats. 

We documented two peaks in the number of megafaunal taxa on 
the HTS, the first of which corresponds to previous studies. The 
interface between warm and cool layers in the Rockall Trough at 
around 1100 m water depth are known to increase benthic abund- 
ance and diversity 7 , both of which peak around 1400 m water 
depth 27 . Likewise, richness on the HTS peaked at 1300-1400 m water 
depth, but we also uncovered an additional distinct peak between 
1500-1600 m water depth. Our CTD data showed transitions from 
WTOW to LSW at these depths, which could foster the rich mixture 



of fauna we observed, some with cool Arctic or warmer Lusitanian 
biogeographic affinities. 

Our ROV and CTD investigations captured a snapshot of the total 
environmental and biological variability on the HTS. Over the last 
two decades, upper and lower layer water mass pH in the Rockall 
Trough has dropped along with the aragonite saturation state 25 while 
levels of DIC have increased, all of which challenge the ability of 
deep-sea species such as cold-water corals to adapt or acclimate 47 . 
The discovery of periodic aragonite undersaturation on the deep 
flanks of the HTS likely has implications for ecosystem functioning 
as well. Organisms with aragonitic calcium carbonate skeletons such 
as S. variabilis are at risk from dissolution and must either be physio- 
logically pre-adapted to these conditions or struggle to produce their 
skeletons. 

Combined with the rapid pace of climate change 25 , deep-sea habi- 
tats in the Rockall Trough are now marked by the footprint of com- 
mercial bottom fisheries, an activity with an impact at least an order 
of magnitude greater than all other human activities in the northeast 
Atlantic 51 . Although we documented 16 sites with elasmobranch egg 
cases, densities of sharks and other fish such as orange roughy on the 
HTS have declined 40-fold since the 1970s 52 and suggest that any 
such spawning grounds should be urgently prioritised as deep-sea 
conservation features. No orange roughy were recorded during our 
survey, the intensive fishery for this species having declined long 
before scientific investigation of this seamount began, as is often 
the case 53 . 
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Knowledge of seamount ecology in the wider northeast Atlantic is 
still poor 54 , but multidisciplinary approaches to collecting these data 
hold much promise in addressing this gap 55 . Our first investigation of 
the turnover in the benthic communities and environmental 
dynamics at the HTS will contribute to the baseline knowledge 
needed to monitor this and other seamounts in the northeast 
Atlantic, while subsequent visual surveys should incorporate detailed 
acoustic data and broader coverage, particularly of the diverse flank 
habitats. 

Methods 

ROV seabed surveys. The first visual survey of the HTS was conducted as part of a 
wider climate change-themed research expedition, the 'Changing Oceans 
Expedition', on the RRS James Cook in May-June 2012 15 . Three dive transects 
(Figure 1, Table 1) were conducted using the Irish Marine Institute's Holland-I, a 
work class Quasar ROV manufactured by Soil Machine Dynamics Ltd. The Holland-I 
was equipped with two manipulator arms, an Insite Mini Zeus video camera with 
direct HDSDI fibre output, a Kongsberg 14-208 digital stills camera (used for the 
present study), two 400 Watt Deep-sea Power & Light SeaArc2 HMI lights, two 
25,000 lument Cathx Ocean APHOS LED lights, and two Deep-sea Power & Light 
lasers 10 cm apart. As there was only time for three ROV transects, dive sites were 
selected to enable preliminary comparisons between the seamount flanks (dives 035 
and 037) summit (dive 036). 

CTD profiling and processing. Oceanographic data obtained by CTD were used to 
help characterise the chemical oceanography of all the different water mass layers at 
the HTS in order to help explain any differences in species richness or composition 
across the seamount. 

A CTD profile was conducted on the seamount along a full-depth cast at 
56°35.295 N and 10°19.65 W. A Sea-Bird 911 plus CTD system (9plus underwater 
unit and Sea-Bird 1 lplus deck unit) fitted with a Rosette water sampling unit and 



10 L Niskin water bottles was used. Conductivity, temperature and pressure 
sensors were laboratory calibrated prior to the cruise, giving coefficients for a 
linear fit. Also attached the Rosette was a Sea-Bird 43 dissolved oxygen sensor, 
Chelsea Aquatracka MKIII fluorometer (set to detect chlorophyll a: excitation 
wavelength of 430 nm and emission wavelength of 685 nm), and Chelsea 
Aquatracka MKIII Transmissometer, to measure dissolved oxygen, chlorophyll 
fluorescence and a particle attenuation coefficient (Cp) measured at wavelength of 
660 nm, respectively, and as described in Findlay et at 23 , Data processing was 
performed using the software SBE Data Processing (V7.21g) and Ocean Data 
View (V4.3.7) for visualisation. Dissolved oxygen (DO) was calibrated against 
Winkler titrations made on water samples taken at various depths through the 
water column as described in Findlay et al. 23 . The offset between titrations and 
sensor measurements was <1%. 

Borosilicate glass bottles with ground glass stoppers (50 mL) were used to collect 
seawater from the Niskin bottles for the dissolved inorganic carbon (DIC) and total 
alkalinity (TA) analysis, at water depths of near surface (2 m), 10, 50, 250, 500, 1000, 
1500 and 1930 m. Sample bottles were rinsed, filled and poisoned with 10 uL mer- 
curic chloride for storage, according to standard procedures detailed in Dickson et 
al. 56 . Duplicate samples were taken from the same Niskin bottle. Samples were 
brought into the chemical laboratory onboard during the cruise at sea, normalised to 
room temperature (approximately 24°C), and analysed for DIC and TA within 
24 hours of collection. 

A Dissolved Inorganic Carbon Analyser (Apollo SciTech, Model AS-C3) calibrated 
using C0 2 Certified Reference Materials (Dickson, Batch 113) was used to measure 
DIC, and duplicate measurements were made to provide an estimate of measurement 
error, which was 0.2%. DIC was also corrected for the addition of mercuric chloride. 

An open-cell potentiometric titration method using an automated titrator (Apollo 
SciTech Alkalinity Titrator Model AS-ALK2) was used to measure TA, calibrated 
with C0 2 Certified Reference Materials (Dickson, Batch 113). Duplicate measure- 
ments were made for each sample to provide a measurement error, which was 0.4%. 
TA was also corrected for the addition of mercuric chloride. 

TA and DIC, together with depth, temperature and salinity, were used to calculate 
the remaining carbonate system parameters (pH, pC0 2 , Q aragonite) using the 
programme C02sys as described in Findlay et al. 23 
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Figure 6 | Species-environment biplots from the flanks (dives 035 and 
037) and summit (dive 036) showing the species that were most strongly 
correlated to the ordination axes. Abbreviations are as follows: Sfrag = 
Syringammina fragilissima; Isid = cf Isidella indeterminate species; Svari = 
Solenosmilia variabilis; Abeat = Aprocallistes beatrxx, Tunic = Tunicata; 
Bris2 = Brisingiid sp. 2; Mling = Mycale lingua; Hdedr = Hexadella 
dedritifera; Grsponge = green encrusting sponge; Hamac = ~<Hamacantha 
indeterminate sp.; Redanem = red anemone sp.; Redsponge = red sponge 
sp.; Redsolcoral = red solitary coral; Orananem = orange anemone sp.; 
Hymed = Hymesdesmia cf (Hymedesmia) curvichela; Aster = asteroid sp. 
1; Redechin = red echinoid sp.; Barnacle = barnacle; Orangesol = orange 
solitary coral; Crinoidl = crinoid sp. 1; Sabyss = Stichopathes cf ahyssicola. 



Best subsets regression analysis using Minitab (v. 16.1.0) showed temperature and 
salinity to be strong predictors for oxygen saturation state using data from the CTD 
profile (r 2 = 99.0, r* M ; = 99.0, Mallows Cp = 48.6, S = 0.32238) according to 
equation (1): 

Oxy %Sat = (3557- 100 * Salinity + 5.83 * Temp) -8.1228; for depth > 1000 m (1) 

Best subset regression analysis of the carbonate system data also revealed that tem- 
perature and oxygen saturation (%) were strong predictors for aragonite saturation 
state (r 2 = 98.0, r 2 ^ = 97.6, Mallows Cp = 4.7, S = 0.082878) according to the 
equation (2): 



^Aragonite = -1.54 + 0.0302 * Oxy%sat + 0.0877 * Temp 



(2) 



Additional CTD data were collected from a Sea-Bird microCAT CTD attached to the 
ROV. CTD data processing were performed using the software SBE Data Processing 
(V7.21g) and Ocean Data View (V4.3.7) for visualisation, and time-matched to the 
latitudinal and longitudinal positions collected from the ROV log. Equations (1) and 
(2) were applied to the ROV CTD data to produce estimates for oxygen saturation (%) 
and 0 Aragonite . 

Image analysis. A subset of ROV still images was created by selecting images taken 
every 60 seconds during the dive. Blurred images, or those taken at different angles of 
the same object, were excluded leaving a reduced subset of 697 images. Each image 
was opened in Picasa viewer v. 3.9 to categorise predominant sedimentary facies and 
record the presence and absence of megafauna. 

For each image, the substratum that was proportionally dominant (>50% of the 
image) was categorised using a modified Wentworth scale 57 (Supplementary Table 
S4). A species matrix was created by enumerating every taxonomically distinct 
organism in an image as present or absent. Given the high uncertainties of identifying 
benthos from ROV images without accompanying physical samples, we adopted best 
practices in habitat analyses 44 to provide the best possible yet conservative identifi- 
cations: taxonomic accounts of the fauna in this and adjacent regions were used to 
distinguish morphologically different species, and the oceanographic setting and 
associated communities were established to support our preliminary identifications 
particularly in habitat -forming fauna such as sponges and corals. Finally, our iden- 
tifications of these and other fauna were verified by senior taxonomic specialists (see 
Acknowledgements). Indeterminate taxa that could not be confidently distinguished 
from species already listed in the matrix were excluded from the analyses. 

Statistical analyses. Trends in species richness were estimated by multivariate 
regression to predict the number of square-root-transformed species (richness) at 
each photographic station using depth, predominant substrata (Supplementary Table 
S4) and oceanographic layer (dynamic WTOW or stable LSW). Richness was not 
normally distributed during any ROV dive (Kolmogorov-Smirnov D dive035 = 0.143, 
p < 0.001; Dji^oae = 0.151, p = 0.039; T> dive037 = 0.135, p < 0.001); thus, non- 
parametric locally weighted scatterplot smoothing (LOWESS) regression was 
conducted for each dive (using XLSTAT v.2008). The standard deviation was used as 
the bandwidth to smooth the data in a second order polynomial function. 

Analysis of similarities (ANOSIM) and similarity of percentages (SIMPER) were 
used to conduct preliminary analyses of environmental variables that correlated to 
changes in species composition. Pairwise measures of the Sorensen Index was used to 
measure community dissimilarity between images, and statistical significance was 
estimated using ANOSIM between categories of geomorphological feature (summit 
or flanks), predominant substrata (Supplementary Table S4), and between oceano- 
graphic layers (dynamic WTOW or stable LSW layer). Each ANOSIM was based on 
999 permutations and yielded an overall global and pairwise values of R that measure 
the strength of the correlation between the environmental variable and species 
turnover. SIMPER was conducted for each category to characterise the megafaunal 
assemblage. 

Following this preliminary investigation, we progressed beyond the correlation 
stage to examine whether environmental variables could explain a statistically sig- 
nificant or large proportion of the species turnover observed during each ROV dive. A 
canonical variance partitioning statistical approach was used to measure the fraction 
of species turnover explained specifically by environmental variability, geographic 
distance between image stations, or induced spatial effects wherein environmental 
dissimilarity varies with the geographic distance between two images. Redundancy 
analysis (RDA) was used to constrain species ordinations into linear combinations of 
explanatory variables, the latter which were then partitioned into the explained by the 
environment, space, spatially structured environmental gradients, and unconstrained 
(unexplained) residual inertia. The constrained variance can be further partitioned 
into the three fractions of pure environmental, spatial and spatially- induced envir- 
onmental effects by conducting partial canonical analyses (pRDA) as described by 
Borcard et at 58 . Response variables (species) were Hellinger -transformed presence- 
absence data 59 . Explanatory variables were depth, the categorical variable oceano- 
graphic layer, and the UTM co-ordinates for each image. Variance partitioning was 
conducted using CANOCO v.5.02 60 . 
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